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Abstract. We perform an analysis of the form factors that rule the structure-dependent amplitude in the 
radiative pion decay. The resonance contributions to n — > ez/ e 7 decays are computed through the proper 
construction of the vector and axial- vector form factors by setting the QCD driven asymptotic properties of 
the three-point Green functions (VVP) and (VAP), and by demanding the smoothing of the form factors 
at high transfer of momentum. A comparison between theoretical and experimental determinations of the 
form factors is also carried out. We also consider and evaluate the role played by a non-standard tensor 
form factor. We conclude that, at present and due to the hadronic incertitudes, the search for New Physics 
in this process is not feasible. 



1 Introduction 



The radiative decay of the pion is a suitable process to be 
analysed within Chiral Perturbation Theory (xPT) [1,2], 
the effective field theory of QCD in the very low-energy 
region. This framework provides the structure of the rel- 
evant form factors through : i) a polynomial expansion in 
momenta, essentially driven by the contributions of heav- 
ier degrees of freedom that have been integrated out, and 
ii) the required chiral logarithms generated by the loop 
expansion and compelled by unitarity. Both contributions 
correspond to the chiral expansion in p 2 /M v and p 2 / A 2 ,, 
respectively, where My is the mass of the lightest vec- 
tor resonance, A x ~ 4irF and F is the decay constant of 
the pion. Hence their magnitude is, in principle, compara- 
ble. The chiral logarithms have thoroughly been studied in 
later years up to 0(p 6 ) both in SU{2) [3] and SU(3) [4,5]. 

However the size of the polynomial contributions is 
more controversial. They involve short-distance dynam- 
ics through the chiral low-energy constants (LECs) of the 
XPT Lagrangian and their determination from QCD is 
a difficult non-perturbative problem. Phenomenology and 
theoretical arguments suggest that the main role is played 
by the physics at the scale My, i.e. the physics of low- 
lying resonances, and a lot of effort has been dedicated 
to pursue this goal [6,7]. This assumption, widely known 
as resonance saturation of the LECs in xPT, implies that 
the structure of the form factors is given by the pole dy- 
namics of resonances and this hint works well in all known 
cases in the allowed range of energies as, e.g. scalar, vec- 
tor and axial- vector form factors in hadronic decays of the 
tau lepton [8,9]. 

Because M^ <C My, it is expected that the structure 
of the form factors in pion decays should be less relevant 
and, accordingly, the approach provided by xPT should 
be good enough even when only the first terms of the 
chiral expansion are included. This is the case of the ra- 



diative pion decay, namely ir — > 1 1^7, I = /i, e, where 
constant form factors (that correspond to the leading or- 
der contribution in the chiral expansion) have widely been 
employed in the analyses of data. However the P1BETA 
collaboration [10] showed that a strong discrepancy be- 
tween theory and experiment arises for the branching ra- 
tio of the process in a specific region of the electron and 
photon energies. Lately the same collaboration, after their 
2004 analysis, concludes that the discrepancy has faded 
away [11]. Curiously enough this decay has a persistent 
story of deceptive comparisons between theory and exper- 
iment [12] that have prompted the publication of propos- 
als beyond the Standard Model (SM) to account for the 
variance [13-17]. Between these it has received particular 
attention the possibility of allowing a tensor contribution 
that could explain the discrepancy by interfering destruc- 
tively with the Standard Model prescription though show- 
ing some inconsistency with the corresponding tensor con- 
tribution in nuclear (3 decay [18]. Related with this issue 
it is essential, seeking to discern the presence of a new 
physics contribution to the radiative decay of the pion, to 
provide an accurate profile of the involved form factors 
within QCD. 

In order to settle the Standard Model description of 
the vector and axial-vector form factors participating in 
7r — * Ivf) decays we study, in this article, the struc- 
ture provided by the lightest meson resonances. This is 
very much relevant on the experimental side because high- 
statistics experiments as PIBETA [11] already are able to 
determine, for instance, the slopes of the form factors in- 
volved in these decays. The procedure consists, essentially, 
in the construction of the relevant three-point Green func- 
tions (GF) in the resonance energy region with a finite 
spectrum of states at the poles of the corresponding mero- 
morphic functions. This can be carried out by employing a 
parametric ansatz or a Lagrangian theory (we work it out 
in both cases). In a second stage several constraints are 
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imposed on the parameters or coupling constants, namely 
chiral symmetry at q 2 <C M R and the asymptotic behav- 
ior ruled by the Operator Product Expansion of the GF 
for q 2 3> Mjj, where Mr are short for typical masses of 
the resonances in the poles and q 2 an indicative squared 
transfer of momenta of this energy region. This method, 
guided by the large-iVc limit of QCD [19] (Nc is short for 
the number of colours), has proven to be efficient [20-26] 
in order to collect the constraints that drive form factors 
of QCD currents. 

Section 2 is devoted to the construction of the vector 
and axial- vector form factors using the procedure outlined 
above. The study of the beyond the SM contribution of a 
tensor current to the radiative pion decay is performed in 
Section 3 and, finally, Section 4 describes the analysis of 
the photon spectrum in this process. After the Conclusions 
in Section 5, three appendices complete the setting of this 
article. 



2 Radiative pion decay : Vector and 
axial-vector form-factors 

The amplitude that describes the ir + — > t + vri process can 
be split into two different contributions : 

m(tt + -> e+vtn) = m ib + m sd ■ (i) 

Here Mjb is the inner bremsstrahlung (IB) amplitude 
where the photon is radiated by the electrically charged 
external legs, either pion or lepton; consequently the in- 
teraction is driven by the axial- vector current. Msd is the 
structure-dependent (SD) contribution where the photon 
is emitted from intermediate states generated by strong in- 
teractions. In this later case both vector and axial-vector 
form factors arise from the hadronization of the QCD cur- 
rents within the Standard Model. 

Because tt + — > e + v e is helicity suppressed, the IB con- 
tribution to its radiative counterpart suffers the same in- 
hibition and, consequently, the electron case is the appro- 
priate channel to uncover the non-perturbative SD ampli- 
tude. Contrarily, the 7r + — > /i + i/^7 decay is fairly dom- 
inated by IB. As a consequence the ir + — ► e + v e ^ is of 
great interest to investigate the hadronization of the cur- 
rents contributing to the SD amplitudes that are driven, 
within the Standard Model, by the vector (Fy(q 2 )) and 
axial- vector (^(g 2 )) form factors defined byQ : 

(7|«7 a d|7r-) = - -r^— £ !i * F v {q 2 ) e af3 ^r» p v , 

(7|u7 Q 7 5 d|7r~) = i ——e^FAiq 2 ) [(r • p)g a& , - p a r ] 

+ ies* a V2F, (2) 

where r and p are the pion and photon momenta, respec- 
tively, q 2 = (r — p) 2 and e is the electric charge of the 
electron. The second term in the matrix element of the 

We use the convention £0123 = +1 for the Levi-Civita ten- 
sor e^vaii throughout this paper. 



axial-vector current corresponds to the pion pole contri- 
bution (which coupling is given by the decay constant of 
the pion F) to the IB amplitude. 

For the full set of expressions for the differential decay 
rate of 7r — > 1 1^7 in terms of the vector and axial- vector 
form factors see, for instance, Ref. [27]. 

Form factors drive the hadronization of QCD currents 
and embed non-perturbative aspects that we still do not 
know how to evaluate from the underlying strong interac- 
tion theory. Their determination is all-important in order 
to disentangle those aspects. It is reasonable to assume, as 
has been common lore in the literature on this topic, that 
hadronic resonance states should dominate the structure 
of form factors and, accordingly, meromorphic functions 
with poles in the relevant resonances coupled to the corre- 
sponding channels have been extensively proposed in order 
to fit hadronic data. This procedure by itself is, however, 
not fully satisfactory because it does not impose known 
QCD constraints. 

On one side chiral symmetry of massless QCD drives 
the very low-energy region of form factors [1]. Hence the 
latter have to satisfy its constraints in this energy domain. 
On the other, one can also demand that form factors in 
the resonance energy region should match short-distance 
QCD properties. This idea was pioneered by Ref. [25] and 
has been used extensively in the last years [20-24,26]. 

In the following we apply these techniques in order to 
determine the Standard Model description of the vector 
and axial-vector form factors in the radiative pion decay. 
Their definition, given by Eq. @, illustrates the fact that 
they follow from three-point Green functions of the cor- 
responding QCD currents. The proper GF in this case, 
namely (VVP) and (VAP), happen to be order parame- 
ters of the spontaneous breaking of chiral symmetry hence 
free of perturbative contributions in the chiral limit. This 
is a key aspect required by our procedure. Hence we con- 
sider in this article their study in the chiral limit, that 
otherwise should provide the dominant features. In the 
following we handle the GF in order to provide a descrip- 
tion constrained by QCD and then we will work out the 
form factors. 

2.1 Vector form factor 

The relevant GF is the (VVP) defined by : 

(nwp)t c (p,q) = (3) 

i 2 J d 4 xd 4 ye^ x+ ^ (0\T{V«(x)V„ b (y)P c (0)} |0) , 
where : 

V«(x) = (V 7 „y V>) (aO , P a {*) = (V^7 5 A» (x) , 

with a = 1,...,8, octets of vector and pseudoscalar cur- 
rents. Moreover SU(3)v symmetry, parity and time re- 
versal demand that : 

(n VVP ); b „ c (p,q) = d abc e^p a q? II VVP (p 2 ,q 2 ,r 2 ), (5) 
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with — (p + q)fj_. In addition Bose symmetry requires 
n VVP {p 2 ,q 2 ,r 2 ) = n VVP (q 2 ,p 2 ,r 2 ). 

The vector form factor Fv(q 2 ), defined by Eq. ([2]), de- 
rives from IIyyp(p 2 , q 2 ,r 2 ) in the chiral limit by : 

Fv(q 2 ) = 2 lim r 2 n VVP (p 2 ,q 2 ,r 2 ), (6) 

with Bo = -(Hj)JF 2 . 

We do not know how to determine 77yyp(p 2 , g 2 , r 2 ) 
(IJvvp for short) from QCD in all energy domains but 
our knowledge of the strong interaction theory allows us 
to know this function in two specific limits as we will com- 
ment now [21,22,28] : 

i) Very low-energy region. The GF has to satisfy the con- 
straints of chiral symmetry encoded in xPT. The lead- 
ing 0(p 4 ) contribution in the chiral expansion is pro- 
vided by the Wess-Zumino anomaly [29] and gives : 

n VVP (p 2 ,q 2 ,r 2 ) ~ * (0 +0(p 2 ,q 2 y^ , (7) 

ii) Asymptotic energy region. The perturbative QCD de- 
termination of the GF, within the Operator Product 
Expansion (OPE) framework gives, in the chiral limit 
and at O{a s )E- 

lim n VVP ((Xp) 2 , (Xq) 2 , (Xp + Xq) 2 ) = (8) 

A— >oo 

_ B F 2 p 2 + q 2 + r 2 fl\ 
A 4 p 2 q 2 r 2 \X e J ' 

lim n VVP {{Xp) 2 ,{q-Xp) 2 ,q 2 ) = (9) 

A^oo 

X 2 p 2 q 2 \X 3 J ' 

lim n VVP ({Xp) 2 ,q 2 ,(q + Xp) 2 ) = (10) 

A — >-oo 

where 
s ab (nvr), pa (p) = 

i J d**e*- (0\T{V«(x)T b pa (0)}\0) , 

( n VT) ppa (p) = (p P 9p,o ~ Po9 w ) n VT (p 2 ) , 

2 Although this result is completely symmetric in the three 
momenta, this does not longer hold when gluon corrections 
are included [30]. Hence we do not expect that this symmetry 
will be sustained when constructing the GF in the resonance 
energy region. 



TpA x ) = (Wy^ (x) , (11) 

where (nvT)iipo-(p) stands for the Vector Tensor GF0. 

In addition to these constraints on the (VVP) GF there 
is also a requirement that we will enforce in any hadronic 
form factor of vector or axial-vector QCD currents. It 
is known [31] that the leading perturbative contribution, 
within QCD, to the spectral functions of both vector and 
axial- vector correlators is constant. Then it comes out, as 
a heuristic deduction, that any of the infinite hadron con- 
tributions to the spectral functions should vanish at high 
transfer of momentum. This implies, in order, that hadron 
form factors of those currents should behave smoothly at 
high energy [32]. Incidentally this feature coincides with 
the known Brodsky-Lepage condition on form factors (de- 
rived within a partonic framework) [33]. Specifically the 
condition, in our case, reads from Eq. ([6]) as : 

lim r 2 n VVP (p 2 ,q 2 ,r 2 ) = . (12) 
p 2 , r 2 — > 
q 2 —>■ oo 

Our task is to construct a function for (VVP) that satis- 
fies, at least, the set of conditions in Eqs. (|7I8I9I10I11I12|) . 

To proceed we use the ideas of the large- Nc limit of 
QCD [19] that, essentially, tell us that GF of currents 
should be described, in the Nc — * oo limit, by meromor- 
phic functions emerging from a theory with an infinite 
hadron spectrum of stable states. This setting is difficult 
to handle and, in practice, one cuts the spectrum reducing 
it to the lightest Goldstone or resonance mesons that obey 
specific hints @. 

In Ref. [22] a Lagrangian theory, including one mul- 
tiplet of vector resonances only, was designed in order to 
obtain an expression for IIvvp that satisfied all conditions 
but for the one in Eq. (TT2"]) . Indeed the fact that only one 
multiplet of vector resonances was not enough in order to 
satisfy all short-distance constraints for this GF was al- 
ready previously noticed [21] with the use of a parametric 
ansatz. It is already well known [21,23] that the MHA is 
more involved if we want that our representation of the 
GF satisfies both OPE and the Brodsky-Lepage require- 
ments. 

The obvious extension is to extend our spectrum by 
including also a multiplet of pseudoscalar resonances in 
the construction of the Uyvp function. Although it can 
be shown that indeed this parameterization satisfies con- 
ditions in Eqs. (|7l8|12p . it fails to meet the OPE condition 
in Eq. ([9]) 0- It is not difficult to relate this problem to the 
fact that the (AP) correlator in the chiral limit is satu- 
rated by one pion exchange [30] . 

It can be seen that all conditions are met if we consider 
in the spectrum of the (VVP) GF two non-degenerated 



3 We use = (i/2) [7 M ,7„]. 

4 This method is known as Minimal Hadronic Ansatz (MHA) 
[20]. 

5 Constraints in Eqs. (|10llip are, in this case, undetermined. 
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multiplets of vector resonances, together with the Gold- 
stone pseudoscalar mesons. Then the ansatz would read : 



n r v ^p(p 2 ,q 2 ,r 2 ) 



(13) 



P + P 1 +P 2 + P 3 



(M 2 



Vi 



q 2 ) (M 2 



Vi 



p 2 )(M 2 2 -q 2 )(M 2 2 ~p 2 )r 2 



with 



Pn = y^^~! Cn-k,k-l,l [r' 2 ) 
k=0 1=0 



2\n-k i „2\k-l 



(9T-V)' 



(14) 



where the coefficients are symmetric under interchange of 
the last two indices Ckmi — Ckim and My 4 i — 1,2 are the 
masses of two nonets of vector resonances in the Nc — > oo 
limit 0. In order to satisfy the short-distance constraints 
several conditions on the djk coefficients arise. Hence the 
chiral symmetry behavior gives : 



cooo = B M^M^ 2T ^ 



The OPE condition in Eq. 



Nc 
4tt 2 

is satisfied for 



c 300 — c 030 — c 120 — C210 

co2i = cm = — Bq F 2 . 







(15) 



(16) 



Finally the Brodsky-Lepage behavior on the vector form 
factor, defined by condition (TT2"|) , fixes one additional pa- 
rameter, namely : 

c 020 = . (17) 

Our ansatz, with all these constraints, satisfies also the 
OPE conditions in Eqs. ([9"TT0"|) . 

If we evaluate now the vector form factor from Eq. ^ 
we get : 



F v (q 2 



cooo + coio Q 



3V2B FM 2 M 2 (M 2 ± - q 2 ) (M 2 2 - q 2 ) ' 



( 18 ) 

and we observe that only one parameter, Coio, has not 
been fixed by our procedure. The expression for the vector 
form factor in the radiative pion decay given by Eq. (|18l) 
is the most general one that satisfies the short-distance 
constraints specified above. As the transferred momenta 
in the 7r + — * e + i> e j process is small by comparison with 
the mass of the lightest vector meson resonance, q 2 <C My, 
it is appropriate to perform the relevant expansion until 
first order in q 2 . Using the result for cooo given by Eq. (fl~5"j| 
it gives : 



F V (q 2 ) =F V (0) 



where Ay 
in 1/Nc and 



1 + A; 



A v 



q 



Mt 



0(q 4 ) 



(19) 



admits an expansion 



Mo) = 



V2N C M^ 
24 Tr 2 ^ 



6 Hadron phenomenology of the pion vector form factor sug- 
gests that the mass of the lightest vector multiplet in this limit 
is very well approximated by the p(770) mass. Hence we will 
take Mvi = M p ^ 70) . 



(20) 



M 2 + M 2 



M 2 Vl 



K 2 



M 2 + £21£ 
cooo 



We must compare our valueQ for Fy(0) ~ 0.0271 with 
the result coming from P(tt° — > 77) and CVC, Fv(0) = 
0.0261(9), and with the recent experimental fit by the PI- 
BETA collaboration, F v (0) = 0.0259(18) [34]. We recall 
that our result for the vector form factor (TT5)) arises from 
a large- Nc procedure where a model of the Nc — > 00 has 
been implemented, namely the cut in the spectrum. At 
q 2 <C My this form factor has been studied up to 0(p 6 ) 
in %PT [3,4]. At 0(p A ) the Wess-Zumino Lagrangian de- 
termines Fy(0) as given in Eq. (|20|) . Higher chiral order 
corrections to this result vanish in the chiral limit, accord- 
ingly their size is suppressed over the leading order by 
powers of M 2 /M v or M 2 /yl 2 that are tiny. Indeed, using 
the 0(p 6 ) odd- intrinsic parity Lagrangian worked out 
in Ref. [35], this modification to i*V(0) is proportional to 
a low-energy constant as M 2 + Cj, that also contributes 
to the 7T° — > 77 decay. From the latter one obtains [36] 



a 



ir 
7 



(0.013 ± 1.17) x 10-^GeV- 



compatible with 



The slope Ay arises at 0(p 6 ) with the usual two 
features : The local operator O22 in provides the 
Nc — > 00 contribution : 



AXPT 



64 7T 2 

~Nc~ 



M 2 + CY 2 r (p), 



(21) 



that is O(l) in the large-Ac expansion, and a one-loop 
calculation that provides the chiral logs corresponding to 
the next-to-leading order [4] : 



A xPT 



M 2 



48 tt 2 F 2 



1 + log 



Ml 

M 2 



(22) 



There is another process directly related with the (VVP) 
GF, namely ir — > 77*; hence it should be related with the 
radiative pion decay. Indeed within the assumptions that 
carried us to Fv(q 2 ) in Eq. (JTSj) , the momenta structure 
for the 7r — ► 77* decay should be the same, though with 
different normalization. In consequence the A]^ ^ slope 
in Eq. (|20)) is the same for both processes. 



The 



7e + e 



amplitude can be expressed by 



77 



>77 



M 1 



(23) 



In the following numerical determinations we will use 
F = 0.0924 GeV, = 0.138 GeV, M K = 0.496 GeV, 



M n o = 0.135 GeV, 
0.775 GeV. 



0.140 GeV and M- 



p(770) 



Vicent Mateu, Jorge Portoles: Form factors in the radiative pion decay 



5 



where q 2 = (p„+ +p e - ) 2 ■ The slope arises at 0(p 6 ) in xPT 
and it is [37] : 



N c 



96tt 2 F 2 



(24) 

Fortunately it has been measured rather accurately [38], 
A 7 = 0.032 ±0.004 and then we can input this measure to 
determine the LEC C^(n) obtaining : 



C% r (M p ) = 7.0±{% x 10^ GeV 



+ 1.0 



(25) 



2.2 Axial-vector form factor 

We now come back to the axial- vector form factor defined 
by Eq. ^j. In order to determine the F A (q 2 ) form factor 
we follow an analogous procedure to the one outlined be- 
fore for the vector form factor. The relevant GF is, in this 
case, the (VAP) defined by : 



(n V Ap) a Jl c (p, q) 



(29) 



i 2 / d 4 xd 4 ye l< - p - x+q - y) (0\T{V^(x)A b u (y)P c (0)}\0) 



where the error includes also the incertitude of the renor- 
malization point fi between M p and 1 GeV. Coming back 
to the slope of Fy(q 2 ) we get : 



Ay = 0.041 



+0.004 
-0.007 ' 



(26) 



that compares well with the recent PIBETA measurement 
A y = 0.070 ±0.058 [34]. 

By comparing now Ajf^^ in Eq. |J5D]) and A^fJ^^ in 
Eq. (|2ip we can provide a determination for the undeter- 
mined parameter in the GF and then give a full prescrip- 
tion for the Fy(q 2 ) form factor in Eq. (TTg)) . For the mass 
of the first multiplct of vector resonances we take M p (j-j ^ 



and for the second M, 



p(1450) 



1.459 GeV 



= (-0.7 ±0.3) GeV" 2 . 

cooo 



(27) 



Notice that the size of this parameter is of the same order 
that the other two terms in 0O . With this result we 

end the construction of the vector form factor in radiative 
pion decays in the large- Nc limit given by Eq. (TT8"]) . 

It is interesting to compare our results with those in 
Refs. [22,39]. As commented above the construction of the 
(VVP) GF in those references was carried out using only 
one multiplct of vector resonances, hence the vector form 
factor in radiative pion decay did not satisfy the constraint 
in Eq. (TT2"|) . With this setting we were able to give a full 
prediction for the leading contribution to the slope Ay, 
namely, 



M 2 

Alf °^°° - My 



4tt 2 F 2 



N C M 2 



(28) 



Using My = M p we got A]^^^ ~ 0.027 to be compared 
with A^ x — 0.028 ± 0.006 from our analysis above. 

The study of the (VVP) GF along the lines outlined 
in this section can also be carried out within a resonance 
Lagrange theory instead of a parametric representation 
as given by Eq. (fT3|) . We collect this procedure in Ap- 
pendix [A] 



8 Notice that one-loop 0(p 6 ) x?T contributions, encoded in 
A^J^ , coincide in 7r — > ev^ and 7r — ► 77* in the SU(3)v 
limit, i.e. for Mk = M v . 



where 



A a Jx) 



i> 7^75 y ^ ) (x) 



(30) 



and the other currents have been defined in Eq. @. The 
structure of this GF is slightly more complicated than 
(VVP) as it involves two scalar functions : 



(n VAP );i c ( P ,q) = 2r bc ( [ B F 2 



(g + 2-g) M <?„ g^ 



q 2 r 2 



+P^(p,q)F(p 2 iq 2 ,r 2 ) (31) 
+Q^(p,q)G(p 2 ,q\r 2 )Y 



with r M = (p + q) p and 



P^(p, q) = inVv - (p ■ q)g l iu 



(32) 



Qixvip, q) = p 2 q^ + q 2 p^p v - (p ■ q)p^q,y - p 2 q 2 g^ 

Then the axial-vector form factor is obtained through 



Bq F p 2 ,r 2 — >0 



lira r T(p , q , r ) 



(33) 



A detailed study of this GF along the line we have per- 
formed above for the (VVP) function was performed 
in Ref. [24]. One of the conclusions achieved was that 
the inclusion of one multiplct of vector, axial-vector and 
pseudoscalar resonances (together with the pseudoscalar 
mesons) was enough to satisfy the matching to the OPE 
expansion of the (VAP) GF at leading order. Moreover 
the analogous to the Brodsky-Lepage condition (fT2|) . in 
this case, was also satisfied, i.e. the resulting axial- vector 
form factor FA(q 2 ) behaves smoothly at high q 2 . We refer 
the reader to that reference for details. Hence we obtain, 
for Nc — > 00 with a cut spectrum : 



F A (q 2 ) 



V2FM^ 



M\ 



My 



(34) 



where M A is the mass of the lightest axial- vector multiplet 
of resonances in the Nc — ► 00 limit. 

At q 2 <C M A we may resort again to xPT [3,5] with 
the expansion : 



F A {q 2 ) = F A (0) 



1 + \a 



M 2 



(35) 



6 



Vicent Mateu, Jorge Portoles: Form factors in the radiative pion decay 



Table 1. Comparison of theoretical and experimental determinations for the low-energy expansion of vector and axial- vector 
form factors. The PIBETA determination assumes that the axial- vector form factor is constant, i.e. it does not consider a slope. 





Experiment [11] 


SU(2) Ref. [3] 


SU(S) Ref. [5] 


Our work 


Mo) 


0.0258(18) 


0.0271 


0.0272 


0.0271 


Xv 


0.070(58) 


0.044 


0.045 


0.041 


F A (0) 


0.0121(18) 


0.0091 


0.0112 


exp. input 


\a 


not measured 


0.0034 


~0 


0.0197(19) 



Both terms, Fa(0) and slope, satisfy an expansion in 
1/Nc, for instance X A = Aj lo ^ 00 + Af /Nc + .... From 
our result above we get : 



F A (0) = V2FM^ 
M 2 + 



M 2 



M\ 



A A 



Mi 



(36) 



Fa{<1 2 ) arises first at 0(p 4 ) with a constant local contri- 
bution from the xPT Lagrangian, namely : 

F^( q 2 ) = 4V2^f(Ll > + Ll Q ). (37) 

The next corrections appear at 0(p e ) in the chiral ex- 
pansion [5]. One of them results from local operators of 
the 0(p 6 ) chiral Lagrangian that, in the chiral limit, only 
contribute to A^i : 



,,(6) | 
M 



A a \N c —>oo 



Mi 



T ' 



a 



78 



2C S 



1 



87 ■ 



-CI 



(38) 



There is also a subleading term, in the large- Nc expan- 
sion, that comes from one-loop diagrams involving the 
0(p 4 ) chiral Lagrangian. However it only affects -Fa(O) 
and it is zero in the chiral limit. The third correction 
is sub-subleading and results from two-loop diagrams 
evaluated with the 0(p 2 ) chiral Lagrangian. The lat- 
ter contributes both to Fa(0) and A^- All local addi- 



tions, F{ i] (q 2 



and A 



(6) 



jv c -»oo, correspond to our result 
in Eq. (f3*B"l) . i.e. Nc — > oo, when LECs are saturated by 
resonance contributions [24]. Though the full 0{p 6 ) chi- 
ral result is rather cumbersome, the authors of Ref. [5] 
have provided a numerical expression for the renormal- 
ization scale /j, = M p . The conclusion is that, in the chi- 
ral limit, subleading contributions to the slope are neg- 
ligible. Notwithstanding it is relevant to emphasize that 
both xPT results of Refs. [3,5] use models to evaluate the 
resonance contributions to the 0(p 6 ) local terms and, ac- 
cordingly, their final conclusion is tamed by this estimate. 

We turn now to give our numerical results. Contrar- 
ily to what happens in the vector case, where the lightest 
vector resonance mass in the Nc — > oo limit is well ap- 
proximated by the /?(770) mass, the axial-vector mass in 
that limit (Ma) differs appreciably from the lightest mul- 
tiplet of these resonances, namely ai(1260). The result 



Ma = V2Mv was obtained in Ref. [25] by imposing sev- 
eral short-distance constraints on the couplings of the res- 
onance Lagrangian. Lately [32] it has been noticed that 
the inclusion of NLO effects in the large- Nc expansion 
points out to Ma < y2My. These results are rather dif- 
ferent from the mass of the lightest axial- vector meson de- 
termined experimentally Ma — 1.230 GeV ~ M a ( 12 6o) [9] 
but it is important to remind that this resonance is rather 
wide. 

Our strategy is the following : we will use the experi- 
mental value of Fa(0) as given in Tab.[T]to determine Ma 
through Eq. (f36|) ; then we provide a prediction for A^ . We 
findET: 



M A = 998 (49) MeV , 
X A = 0.0197(19) , 



(39) 



where the error stems only from the experimental uncer- 
tainty in Fa(0). Notice that this result satisfies Ma < 
V2M V ^ 1096 MeV. 



2.3 Theory versus Experiment 

We are now ready to compare our results with other 
theoretical settings and experimental determinations. In 
Tab. [T] we compare our outcome for the low-energy expan- 
sion of the form factors with the one provided by 0{p e ) 
xPT and the recent PIBETA published values. 

As Fy (0) is ruled by the Wess-Zumino anomaly all the 
theoretical results agree for this parameter. Leading cor- 
rections to this value are driven by the pion mass and as 
a result happen to be tiny [22]. This is also reflected in 
the excellent comparison with the experimental determi- 
nation. The agreement is also good for the slope of the 
vector form factor, considering the large error of the ex- 
perimental value. 

The axial-vector form factor does not arise a similar 
consensus. As indicated above %PT can only predict re- 
liably all loop contributions (up to 0(p 4 ) in the even- 
intrinsic parity and 0(p 6 ) in the odd-intrinsic parity sec- 
tors) while higher order loops involve the couplings of local 
operators. Moreover tree-level 0(p 4 ) (|37| and 0(p 6 ) (|38l) 



9 It is important to notice that the value of -Fa(O) measured 
by the PIBETA experiment assumes no slope for the axial- 
vector form factor. We should repeat this exercise when Aa is 
included. 
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terms can only be determined in different models for res- 
onance saturation contributions. The excellent agreement 
between the xPT results and the experimental determi- 
nation of Fa(0) is, indeed, not a major issue as the axial- 
vector form factor in radiative pion decay is the main phe- 
nomenological sourc^E] to fix the value of L\ . It happens 
that F%\o) arises from a strong cancellation between the 
Lg and L\ LECs and, in consequence, it is very sensi- 
tive to the chosen value for L\q. In terms of resonance 
saturation this sensitivity moves to the value of the axial- 
vector mass Ma input in the numerical determination. 
The value of L\ ~ —5.5 x 1CP 3 , used by Ref. [5], arises 
for M A ^ 1 GeV. 

Our model of large- Nc gives the leading result for the 
axial-vector form factor parameters and there are lead- 
ing Goldstone-mass driven contributions that we have not 
considered. In the x?T framework these 0(p 6 ) corrections 
arise from the LECs and, a priori, it is difficult to estimate 
their contribution due to our lack of reliable knowledge on 
those low-energy couplings. However it has been pointed 
out [5] that the role of LECs is unimportant in the 0(p 6 ) 
corrections. As the subleading loop contributions are also 
tiny, it is concluded that is not sizeable and Fa(0) is 
ruled by the leading 0(p 4 ) contribution by far. 

However using as input the experimental value of 
Fa{0) we find a large value for Xa- As subleading 1/Nc 
loop contributions seem to be tiny our leading result shows 
a clear discrepancy with the estimates of tree-level contri- 
butions performed in the chiral framework [3,5]. It would 
be very much interesting to have an experimental determi- 
nation of A,4 in order to disentangle the different resonance 
models. 



3 Beyond SM : Tensor form factor 

As pointed out in the Introduction, the history of the ra- 
diative decay of the pion accumulates a few clashes be- 
tween Theory and Experiment. It seems though that, after 
the latest analysis by the PIBETA collaboration, the land- 
scape has very much soothed. However it has become cus- 
tomary to investigate possible contributions beyond the 
Standard Model in order to appease alleged discrepancies. 
Between the latter the possible role played by a tensor 
form factor has thoroughly been studied [13-17]. 

The new short-distance interaction can be written in 
terms of quark and lepton currents and it reads : 

C T = ^=V ud F T [qa IM/ (l- 75 )q] [la^ (1 - 7 s) v t ] , 

~_ X " (40) 
where Ft is an adimensional parameter measuring the 
strength of the new interaction. As the product cr^ 7 5 
is not an independent Dirac matrix (due to the identity 



1 Lg is rather well determined from the phenomenology 
(squared charge radius of the pion) and its numerical value 
agrees nicely with resonance saturation. 



a ^ 7 5 = _ « E «S CTq/3 ) we can wr it e (gO]) as : 

C T = -^V ud F T [qa^ vl5 q] [£a^ (1 - 75 ) v t ] • 

(41) 

In the Standard Model the later structure, a tensor-like 
quark-lepton interaction, arises from loop corrections to 
the tree-level amplitudes and gives a tiny value for Ft ~ 
10 -8 [14]. More sizeable contributions could come from 
New Physics models. Leptoquark exchanges, for instance, 
could give Ft ~ 1(U 3 [15], while SUSY contributions pro- 
vide Ft ~ 10~ 4 — 10~ 5 [14] for light supersymmetric part- 
ners. 

The hadronization of the tensor current, at very low 
transfer of momenta, is driven by the constant fx defined 
by: 

(l\ua^j 5 d\TT-) = --faiPntv-PuCn), (42) 

where p is the photon momentum 0. The determination 
of /t involves QCD in its non-perturbative regime and, 
consequently, is a non-trivial task. We will come back to 
this issue in the next Subsection. 

It is possible to obtain the product T = Ft /t from the 
analyses of different processes. Hence from some previous 
discrepancy in the 7r — > e v e ^ process it is found that 7^ = 
— (5.6 ± 1.7) x 10 -3 [13], while from the introduction of a 
Gamow- Teller term in the amplitude of nuclear /3-decay 
[18] gives T N = (1.8 ± 1.7) x 1CT 3 . 

3.1 (VT) Green function : The tensor form factor 

If we want to extract information on the value of Ft from 
experimental data, we need a reliable QCD-based deter- 
mination of the hadronic tensor form factor. Using LSZ 
and at leading order in the pion mass we can express the 
matrix element (|42|) as follows E3 : 

(7 lua^u 75 d\ 7r~) = -= (7 \wfj, u u + da^d\ 0) = 
y2F 

y/2e 

- n VT{^){p^ -Pvty), (43) 

where in the last step we have used again the LSZ re- 
duction formula applied to the (VT) correlator defined in 
Eq. (JTTJ> . Then we have : 

f T = i^y n VT (0) . (44) 

To determine the (VT) correlator it suffices to employ 
Eq. (p~0|) with our ansatz in Eq. (| 13|) . It gives : 

rr 1 2\ i c no + c 200 + 2c m g 2 

nVT{q) = 2 (M^H^ Va -^) ' (45) 

11 There is in fact another Lorentz structure contributing to 
this matrix element but it carries higher orders in momenta. If 
the latter is included /t acquires a dependence in the squared 
of the transferred momenta, i.e. fr(q 2 )- See Appendix iBl for a 
detailed evaluation of both form factors. 

12 See Appendix[C]for a derivation of this expression. 
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that satisfies the proper high-energy behaviour ruled by 
Eq. (fT0|) . Unfortunately Ilvx(0) is not fully specified by 
short-distance constraints within our approach as we have 
only fixed the cm parameter. 

In order to provide an estimate for fx we consider 
this correlator including one multiplet of vector resonances 
only. To proceed we can use the results for the (VVP) GF 
as in Ref. [22] (using Eq. (jTOJ) for instance) or, equivalently, 
from the (VAP) GF [21,24]. Both procedures yield the 
same result, namely : 



n VT (q 2 



that matches the OPE result [42] 

. B F 2 



2„2\ 



lim IIyx{\ q 

A — >OG 



o 



A 4 



(46) 



(47) 



when the large momentum limit is taken. Notice that the 
result (|46|) can be recovered from (|45f by taking the limit 
My 2 — > co, demanding cm = and identifying the rest 
of the constants. Using Eq. (jl6"| in Eq. (14"4")) yields : 



h 



2y/2B F 
3M 2 



(48) 



An educated guess can be obtained by writing 
B F = - (Hj) /F and the estimate ($ip) (1 GeV) = 

- (242 ± 15MeV) 3 [43]. We obtain fx = 0.24 ± 0.04. 

Another parameter of interest is the susceptibility of 
the quark condensate \ z defined by the vacuum expec- 
tation value of the tensor current in the presence of an 
external source Z^ v [44, 45] : 

(0|^V|0} Z = g+x* {H>) -V- (49) 

In our case we consider the magnetic susceptibility \ given 
by an external electromagnetic field as : 

( 7 I ua^u + da^d | ) = -ie (e„ + e d ) \ (^V , ) F^ v 

with e u — 2/3 and = —1/3. Using the first equality of 
Eq. dm we get : 



(51) 



h = - ^yXB Q F, 

and comparing with Eq. (|48[) we obtain 
2 



X 



M 2 



-3.3 GeV" 



(52) 



There are several determinations of the magnetic suscep- 
tibility that provide a range that runs from x = —(8.16 ± 
0.41) GeV -2 [44] up to x ^ -2.7 GeV" 2 [46]. 

3.2 Lattice data and sum rules 

The last years have witnessed and increasing attention 
to the determination of matrix elements of tensor quark 



currents. For instance, together with the QCD sum rules 
technique [46,47], lattice has also performed evaluations 
of amplitudes involving the tensor current and a vector 
resonance [48,49] : 



(0\V«\p b (p,\)) 



1 



5 ab M v fv< 



V2 



(53) 



(0 | T; v I p\ Pl A)) = --j= S ab (e> 



where p b (p,X) is a vector resonance with momentum p, 
hclicity A and polarization vector e^. The vector and ten- 
sor quark currents have been defined in Eqs. (|4llip . The 
scale dependence of the fy{p) in Eq. (fS"3"]) reflects the fact 
that the tensor current has a non-vanishing anomalous 
dimension. 

Within the large- Nc framework it can be shown [42,50] 
that if we consider a single multiplet of vector mesons we 
get : 



n VT (q)--2 M 2_ q2 ■ 

Comparing with Eq. (|46p we obtain the relation : 

2Bn F 2 



fv fv — 



Ah 



(54) 



(55) 



The fv coupling can be obtained from the measured 
r(p° -> e + e~) [381. We obtain f v ~ 221 MeV with an 
expected tiny error Then from Eq. (|55[) and using the 
value of the quark condensate quoted above we get : 



/^(lGeV) = 165 ±31 MeV, 



(56) 



where the error collects only the uncertainty in the value 
of ^^)_. Our result is in excellent agreement with those 

coming from QCD sum rules : f^- = 160(10) MeV [46] and 
f^ = 157(5) MeV [47]. 

Lattice evaluations determine the ratio with the vector 
coupling. From our results we get 



^(lGeV) 
Jv 



0.75 ±0.14. 



(57) 



to be compared with the quenched value [48,49], run down 
to n = 1 GeV : 



-^-(lGeV) = 0.74 ±0.03. 
Jp 



(58) 



Finally from the later result and the phenomenological 
value of fv lattice provides the determination : 



/^(lGeV) = 164 ±7 MeV, 



to compare with our figure in Eq. 



(59) 



13 The vector coupling can also be determined from short- 
distance analyses within resonance theory [24], giving f v — 

' ' W V, that translates into fv = 207 (15) MeV for Ma = 



Mi —M 



998 (49) MeV (|39[) , in excellent agreement with the quoted phe- 
nomenological result. 
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Table 2. Comparison of the theoretical predictions and the experimental data for TZq = 10 Rq for constant form-factors and 
different predictions of the q 2 dependence. 



Efj n (MeV) 


E™ in (MeV) 


amiii 
e -y 


ft cxp [34] 


IZth (without slopes) 


^th (with slopes) 


1Z th SU(2) [3] 


TZ th SU(3) [5] 


50 


50 




2.614(21) 


2.78(38) 


2.81(38) 


2.46(35) 


2.72(38) 


10 


50 


40° 


14.46(22) 


14.81(54) 


15.08(58) 


14.73(53) 


15.00(57) 


50 


10 


40° 


37.69(46) 


38.08(98) 


38.41(103) 


37.51(94) 


38.17(103) 



4 Analysis of the photon spectrum in the 
radiative pion decay 

The PIBETA experiment has thoroughly measured the 
photon spectrum in the radiative decay of the pion [10]. 
Though the results of that reference seemed to confirm 
a serious discrepancy with theoretical determinations, an 
ensuing analysis of more data and the refinement of sys- 
tematic errors [11,34] has brought a close agreement be- 
tween Theory and Experiment. 

The experimental available data amounts to the 
branching ratio of the radiative pion decay integrated in 
different subregions (Q) of the final state phase space : 

Rq = j^— I dQ 3 £ \M(E e , E v )\ 2 , (60) 

where the sum runs over the polarizations of the final par- 
ticles. The three regions and the experimental results are 
shown in Tab. O 

We test the predictions ruled by our determination 
for the hadronic form factors with the experimental data, 
ignoring first a possible tensor interaction, and compare 
them with other theoretical settings. In order to achieve 
the accuracy required by the experimental information 
higher order radiative corrections to the decay [51] must 
be included and they have been implemented in our anal- 
ysis. The numerical input for vector and axial-vector form 
factors is given in Tab. [1] 

In the fourth column of Tab. [2] the latest experimental 
data are given; in the fifth and sixth we show the results 
provided by our analysis. We study the numerical impact 
of the momenta dependence of the form factors by setting 
the slopes to zero and we conclude that it is tiny : the 
q 2 dependence tends to increase the central value of R 
but the modification is by far within the errors. The last 
two columns bring the results yielded by two- and three- 
flavour two-loop %PT calculations. The evaluation of the 
errors for the theoretical predictions is ruled by those in 
the form factors. The estimate of the latter has been done 
in the following way : We assume no error coming from 
the slopes (since their numerical impact is very poor); to 
the vector form factor we assign the same error as that 
of the experimental determination ~7% and to the axial- 
vector form factor we attach the error of the experimental 
input. Finally the error given for the xPT calculations 
only considers the scale dependence that is a tiny 5%. 

We conclude then that the corrections induced by the 
q 2 dependence of vector and axial- vector form factors are 



numerically negligible unless the theoretical error is re- 
duced. For this we would need a better determination of 
vector and, specially, axial-vector form factors at q 2 = 0. 
When comparing our results with experimental data, we 
see that our predictions are in agreement with previous 
estimates. 

As a final exercise we use the experimental data to fit 
the value of T = FtJt defined above. In order to reach 
this purpose we use the experimentally fitted values for 
the hadronic inputs iV(0) and -Fa(O), and our results for 
the slopes Ay and A a- Finally, to extract the value of the 
Ft coupling from the fit, we use our determination for the 
tensor form factor fx- The value that we obtain is com- 
patible with zero and its order of magnitude is compatible 
with that dictated by SUSY : 

F T = (1 ± 14) x 10~ 4 . (61) 

5 Conclusions 

Radiative pion decay has been a continuous source of de- 
bate between theoretical predictions and experimental de- 
terminations. Nevertheless the latest analysis by the PI- 
BETA Collaboration seems to bring a close agreement be- 
tween both sides. 

In this article we have performed a detailed analy- 
sis of the structure-dependent amplitudes contributing to 
7r — ► ev e ^. The q 2 dependence of vector and axial-vector 
form factors, driven by the Standard Model, has been 
rigorously constructed through the study of the (VVP) 
and (VAP) Green functions, by matching meromorphic 
ansatze with their leading OPE contributions. Moreover 
we have also required that our form factors are soft at high 
transfer of momenta. Hence we obtain the most general 
(and simple) functions that satisfy all those constraints. 
The appropriate structure of the form factors requires a 
double vector resonance pole for the vector form factor and 
a single axial-vector resonance pole for the axial-vector 
form factor. After a small momenta expansion we com- 
pare our results with those of xPT and while in the vector 
sector we find complete agreement, our slope for the axial- 
vector form factor is much larger than the one provided 
by modelizations of local terms in the chiral framework. 

The role of a tensor contribution to the radiative pion 
decay has customarily been taken into account in order to 
analyse the experimental results. We use those in order to 
fix the size of the contribution and we find that it is com- 
patible with zero. Incidentally we have given a prediction 
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for fy that measures the coupling of a vector resonance 
J PC = 1~~ to the tensor current. Our results agree well 
with determinations from QCD sum rules and quenched 
lattice. 

We conclude that the Standard Model is able to em- 
body the experimentally known features of the radiative 
pion decay. As it happens with other decays involving non- 
perturbative strong effects, the rather large size of the nu- 
merical incertitudes generated by our lack of knowledge 
of this QCD regime shows that this process is, at present, 
unsuitable for the search of New Physics. 
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A (VVP) from a Lagrangian 

We used the meromorphic ansatz in Eq. (| 13(1 in order to 
determine the vector form factor in radiative pion decay. 
To reach the same result one could also proceed starting 
with a Lagrangian like the one given by Resonance Chiral 
Theory (R%T). This study, with one multiplet of vector 
resonances, was already carried out in Ref. [22]. However 
we have concluded that QCD constraints seem to indi- 
cate the need of including a second vector multiplet. In 
this appendix we construct the short-distance constrained 
(VVP) function within R%T. 

In order to proceed we need to build the odd-intrinsic 
parity RxT Lagrangian with two multiplets of vector res- 
onances. The interaction pieces are : 

Fv (vrf +lu .) + P^=(vru 



£cvcn 
V 



£odd 
V 



2\/2 

£-wz + i {< 



2V2 

w 



(AT) 



X-f+ /- 



+ 



-w 



£ 



+ iC 22 
7 



v A /£*{tfV})} 



Mv: Vijp 



£ 



My 2 UV * JP 



4 
i=l 



V1V1P 



v 2 v 2 p 



E d nOy lV2 p +dfO Vl y 2 j, 
n—a,b,c,d,e 



(A.2) 



where £\yz is the Wess-Zumino Lagrangian [29] that 
arises at 0(p 4 ) in x?T. As specified above only two opera- 
tors contribute at 0(p 6 ) [35]. Operators Oy i jp and OviViP 



were already given in Ref. [22] and will not be repeated 
here. For the last part of the Lagrangian there are two 
subsets of pieces [40] : 

— V1V2P terms, which contain vertices with Goldstone 
and two vector resonances from different multiplets : 



u ViV 2 P 


— E-fivpa 


({vr,v 2 nv a 


u°) , 




u ViV 2 P 


— E-fivpa 


avr,vr}v a 


«"> , 




^V X V 2 P 


— Efivpcr 


av a vr,v 2 pa } 


u") , 




u VxV 2 P 


— E-fivpa 


({v a vr,vr} 


«"> . 






— £p,i>pcr 


({v°vr,v 2 pa } 


u a ) . 


(A.3) 



V1V2J terms, with two vector resonances from different 
multiplets and one pseudoscalar external source : 



O 



f 

VxV 2 J 



iep Upa ({Vr,V 2 p °} X -) 



(A.4) 



The result for the (VVP) Green function defined in 
Eq. © is : 



tjRxT _ R 
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where 



A' = 8F^ 2 (d[-4) 
B = 8F*d 3 , 



A = 8F V (di-tfe) , 

C = 64F v d 2 , 
C 1 = MF v 2 d' 2l 



D = - 
D' = — 



V "2 ■ 
32V2 F V C 3 

M Vl 
32V2 F v c' 3 



M V2 



4V2F V 
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4 ^(4 + 4-4) 



E 
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M V2 
4^2 F v 



M Vl 
4V2F, 

M V2 
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(A.6) 



In terms of the parameters of our ansatz in Eq. (|13p we 
obtain : 
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My„M- 



Vi 



4tt 2 ' 



(A.7) 



in units of — Bq. Chiral symmetry, implemented in our 
Lagrangian, brings features that with the ansatz had to 
be forced by hand. In this way we immediately find that 
C300 = an d C210 = 0. Moreover, as a bonus we also find 
C200 = 0. 

The rest of constraints are given in Eqs. (|16I17[) . In 
addition we find five more relations : 

w 



D 



D' 
2{K 



G 7 


-0, 


— w 




G22 


= 0, 


G + G' 


= 0, 


E + E' 


-0, 



K>) = -4^ 



4tt 2 



(A.8) 



After applying all the constraints coming from the OPE 
expansion and Brodsky-Lepage asymptotic condition we 
obtain the following relations among the Lagrangian cou- 
plings : 



4c 3 
4 4 



Cl 

4 



Cl 

3t4 
C2 + C 5 

4 + 4 



0, 
0, 
0, 
0, 



F v M Vl M Vl N c 

c 5 - C 6 + ~ ~ (C 5 - Cg) 



F V M V2 

8F 2 d 3 + 8F{?4 + 8F v F v (d c - d d ) 

r- , My - My , 

-8V2F V v \„ ^ (4 5 -c 6 ) 



F v 64V2tt 2 



Mi 



V-2 



+MV ^- F ' 



AF v (di +8d 2 )+ 4F( 2 (d' 1 + 84) 
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-4Fy F v (d a + db — dd + 8d f) transforming with the compensating field 

i 

r * 

M V2 



2 2 To the lowest order in the chiral expansion we have one 

+-^Vi iQ n 2 = ^ ' operator that contributes to both the radiative pion decay 



(A.9) 



and the (VT) function, 



8M V2 F v d 3 + 8F^M Vi d^ £ xPT ± Ai ^ f+ ^ (R9) 



+4F V F V [M Vl (d b + d c -d a - 2d d 



N c 



When explicitly including resonances the corresponding 



+M V2 (d a + d c — d b )] = -M Vi M V2 -^-^ . interacting Lagrangian of interest reads 



+V2F VT M v (V^tt?) . (B.10) 



B Chiral Lagrangians with external tensor 
sources 

The couplings Fy and Fyx are related with those defined 

In this appendix we determine the matrix element in m Eq. ([53]) by f v = V2F V and fy = V2F VT . Upon 

Eq. g2D by incorporating tensor sources to the resonance integration of the vector meson we can relate the couplings 

Lagrangian. The inclusion of those external fields in x?T Q f the two Lagrangians : 
has been studied in Ref. [41]. 

We extend the QCD Lagrangian to accommodate an A —A~ FyFyr y^. 

external tensor source : l — l • \ ■ ) 

Cqcd — £qcd + £ext 5 (B.l) If resonance saturation of the chiral LECs would hold 

- T/ii us i ?/ • si t -u v i when including external tensor sources, yli should be given 

C ext = V7„K + Tfeo")^ " ^ ~ *75P) V + i> • only by the second term in Eq . {HUH). In fact this is 

„. . . . , , . , . , , the case when we enforce short-distance constraints. If wc 

1 he tensor source includes both octet and singlet currents , ,, /TArn \ ^ r , n , • ^ tttk 

evaluate the (VI) Green function defined m Eq. yjj we 



t""=Ef^ (B - 2) 



a=0 

One finds that 



obtain 



M v -p 2 



that satisfies the OPE constraint in Eq. (fTTj) provided 
^ ^ ^ = <Al cr"" V>fl + &i V , (B.3) that : 



and the change of basis reads 



Ai = 0, 

FvFvt = ^^. (B.13) 



IfiV _ r>puXp . , r,PvXp A - i- - i- J , , 

c — i-Ap -t- t Ap , m 

t^ v — p^ x P t Xp , (B.4) Using now the relation in Eq. (|44|) we get again the result 

in Eq. ggj) for / T . 

where P^ v ^ p are the chiral projectors for the tensor fields, As mentioned in Section [3] there are two form factors 
given by involved in the hadronic tensor matrix element 

pp,vXp _ \i p,X up _ uX p,p _ : p,vXp\ <7l"^75d|0= ~ \ frtf) (p p^v ~ Pu^p) (B.14) 

e :g T (q 2 )[e-q (p»q v - p u q») + q-p (q»f? - , 



ppuXp = (pj^M _ (B 5) 



where p is the photon momentum and q the transferred 
In order to manta in th e chiral invariance of the extended momentum by the tensor current. To generate the second 
QCD Lagrangian (|B.1|) the tensor source must transform Lorentz structure one needs operators of higher order in 

the chiral expansion, such as 



tp,v -> gRtp, v g L - (B.6) 



^93 = (V„tf V a / +cw ) , (B.15) 



It is convenient, in order to build an effective Lagrangian 

invariant under the chiral group, to define the tensor op- f Ref. [41]. In the framework of resonan ce chi ral theory 

erators : one needs in addition of the operators in (|B.10[) the basis 

= w t^ v \v ± ut^'u, (B-7) of odd- intrinsic parity operators of Ref. [22]. These new 
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contributions do not modify result of the /x = /t(0). The 
result reads : 



f T (q 2 ) 



9r{q 2 ) 



VT 



3 F M- 



2F V - (4V2My (ci - c 2 - c 5 + 2c 6 ) 

„2 



8 Fyx 

3F My 



+8F V d 3 



2 My 



My 



(B.16) 



My- q ^'^ 



My 



My - q 2 
+V2 



c 6 - 2 c 7 



(B.17) 



F v ( Ml 



M v \M V - q 2 



The spectral function of the tensor-tensor currents cor- 
relator to which the amplitude in Eq. (|B.14[) contributes 
behaves as a constant at high q 2 and leading order in 
as [42]. Hence the fr{q 2 ) and gr(q 2 ) form factor should 
exhibit a smooth behaviour, vanishing at large transferred 
momentum. 

Imposing that the /t(<? 2 ) form factor vanishes at large 
momentum we get the constraint : 



ci - c 2 - c 5 + 2 c 6 



V2 F v 



4 M v 

The same procedure with gr{q 2 ) gives : 
c\ — c 2 — C5 — 2 c-j = 



(l + 4d 3 ). (B.18) 



M v 



(B.19) 



Interestingly enough these constraints fully determine 
both form factors : 



Mq 2 ) 
9r{q 2 ) 



2V2 F VT F V M v 



3F My-q 2 ' 
h(q 2 ) 

My ■ 



(B.20) 



Notice that the contribution of the gr(q 2 ) form factor to 
the matrix element under study is fairly suppressed, typ- 
ically 0{q 2 /My) over the /t(<7 2 ) contribution. 



C LSZ formula for a soft pion 

In this appendix we discuss the derivation of Eq. (|4"3"|) . Let 
us start defining three quark currents : 

A»{x) = d(x)^j 5 u(x), 

Tr(x)=u(x)a»» l5 d(x), (C.l) 
T%" = u(x) a** u{x) + d(x) a» v d{x) . 

Then we can construct the following GF that after partial 
integration can be expressed as : 



(j(p)\d a A*Tr\0) 



(C.2) 



l J d 4 xe-( 7 b)|T{9 Q ^(x),rr(0)}|0) = 
r a [ dSe^< 7 (p)|T{^(aO,Tr(0)}|0) 



i (7(p)|Tf(0)|0) . 



(C.3) 



We can relate this GF to the (^(p) \u ct m „ 75 d\ n (r)) ma- 
trix element (14TJ1) through the LSZ formula : 

(j(p)\ua^ l5 d\7r-(r)) = 
r 2 - M 2 



lim 



f- (j(p)\d a A a Tr\0) ■ (C.4) 



r 2^ M 2 ^/2FM 2 

Then the (j(p) \d a A a Tg v \ 0) GF must have a pole at r 2 = 
M 2 . Since the second piece in (|C.3j) has no r-dependence 
it cannot have a pole, and this must come from the first 
term : 



r a I d 4 xe-(7(p)|r{^( a; ),rr(0)}|0) 



^(p,r). 



1 



r 2 _ M 2 



(C.5) 



As the divergence of the axial current is zero in the chiral 
limit, the left hand side must vanish for M„ — 0, what 
implies : 

^"(P, r)U.=o = * r2 (-Tip) \ T rm 0)m„=o • (C6) 

Finally, assuming that pion mass corrections are small, 
(IC.4[) yields the desired result : 

(j(p)\ua^ 75 d\-K-(r)) = i-^ ( 7 (p)|Tr(0)|0) , 

(C.7) 

that is also the result that stems from the soft pion theo- 
rem. 
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